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Rice  University 
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I*  Introduction 

An  array  of  acoustic  sensors  is  placed  in  &  known  spatial  pattern  to  record  the 
acoustic  environment*  Measurement  of  the  acoustic  field  by  an  array  offers  two  basic 
improvements  over  the  $ignal~processing  capabilities  of  a  single  sensor*  The  first  is  the 
determination  of  the  bearing  of  the  acoustic  source(s)*  Bearing  cannot  be  obtained  with  a 
single  sensor  whereas  an  array  offers  some  bearing  resolving  capability*  This  capability 
is  usually  measured  by  the  just-de tec  table  separation  of  two  equi-strength  sources  for  a 
given  signal-to-noise  ratio  at  one  sensor  in  the  array[6]*  The  determination  of  the 
bearing  of  a  remote  source  of  acoustic  energy  remains  one  of  the  fundamental  problems  of 
passive  sonar  systems*  The  waveforms  recorded  by  each  sensor  can  be  acquired  with  as  much 
fidelity  as  desired*  In  contrast*  the  size  and  geometry  of  the  array  is  usually  limited 
by  physical  considerations;  these  limitations  restrict  the  spatial  resolution  of  source 
bearing*  The  second  improvement  is  the  increase  of  signal-to-noi se  ratio*  If  the  noise 
field  is  oncorrelated  at  each  sensor  location  with  respect  to  all  other  locations*  the 
signal~to-noise  ratio  in  the  array  output  is  increased  by  a  factor  equal  to  the  number  of 
sensors  comprising  the  array*  This  factor  decreases  when  the  noise  field  is  correlated  at 
the  sensor  locations*  The  measure  of  improvement  in  the  signal-to-noi se  ratio  is  array 
gain:  the  ratio  of  the  signal-to-noise  ratio  at  the  array  output  to  that  obtained  with  a 
single  sensor* 

Adaptive  beamforming  methods  (ABF)  are  known  to  have  superior  bearing  resolution 
capabilities  when  compared  to  conventional  beamform  ers  (in  a  theoretical  sense)[61* 
Specifically*  the  minimum  energy  method  has  been  analyzed  extensively  in  this  regard. 
However*  there  is  no  theoretical  basis  indicating  that  this  method  has  the  best~possible 
bearing  resolution  properties*  On  the  other  hand*  the  array  gain  provided  by  this  method 
is  optimum:  no  other  beamforming  technique  can  yield  a  larger  increase  in  signal-to-noi se 
ratio*  This  report  is  concerned  with  a  new  ABF  scheme  which  is  similar  in  many  respects 
to  the  minimum  energy  method*  It  can  demonstrate  greatly  increased  bearing  resolution 
properties  but  at  the  expense  of  array  gain*  The  scheme  is  based  on  an  eigenvector* 
eigenvalue  decomposition  of  the  empirical  correlation  matrix*  which  is  then  truncated  so 
as  to  retain  only  those  terms  which  best  contribute  to  increased  bearing  resolution*  This 
approach  is  similar  to  those  described  by  Schmidt[13]  in  his  MUSIC  system*  by  0wsley[12] 
in  his  modal  decomposition  approach*  and  by  Bienvenu[2*  3]*  Analytic  results  are 
presented  here  which  contrast  the  bearing  resolving  properties  of  these  various 
eigenvector  methods  and  the  minimum  energy  method* 

II*  Preliminaries 

L*t  xa(t)  denote  t he  output*  takes  from  as  array  of  seasors  having  a  kaowa  geometry . 
ooaeists  of  eoapatiag  the  quantity 
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y(t) 


a  x  (t-x  ) 

a  ■  a 


(1) 


where  a  is  the  amplitude  weighting  (shading)  applied  to  the  a^-sensor  ontpnt.  x  is  the 

delay  applied  to  the  tensor  ontpnt#  and  M  is  the  number  of  sensors  in  the  array.  The 

parameters  {a  }  and  {r  }  of  a  beam  are  chosen  according  to  some  desired  criterion  (e.g.# 

m  m 

steering  the  beaa  in  a  particular  direction,  ainiaizing  sidelobe  height#  etc.). 

Beaafoxaing  can  also  be  viewed  as  a  type  of  aultidiaensional  spectral  analysis(l,  9]. 
Evaluating  the  Fourier  transform  of  the  beaa  y(t)#  we  have 

fcl  -j2nfx 

Y(f)  »  )  a  e  "x  (f).  (2) 

JsQ  “  ■ 

Therefore,  at  each  tea^oral  frequency  f,  the  Fourier  transform  of  the  beaa  can  be  written 
as  the  dot  product  of  two  vectors 


Y(f)  -  A'X 


(3) 


where  A  denotes  the  steering  vector  consisting  of  the  elements 

+j2nfx 

a  ® 

A  *  a  e 
a  a 


and  J  denotes  the  vector  comprised  of  the  Fourier  transforms  of  the  sensor  outputs*  Here, 
2'  denotes  the  conjugate  transpose  of  2*  Assume  that  we  have  a  linear  array  of  equally- 
spaced  sensors;  in  this  instance#  the  delays  will  be  of  the  fora  tffl  *  aT.  Equation 

(2)  becomes 


Y(f,  T) 


I 

■H) 


(4) 


Consequently,  Y(f,  T)  is  the  Fourier  transform  of  the  sequence  a  X  (f).  For  a  linear 

a  a 

array,  computing  a  transfoxm  along  the  index  m  is  identical  to  coa^uting  the  transform  in 

space  across  the  array.  Y(f,  T)  is,  therefore,  the  spatial  transfoxm  of  a  X  (f)  evaluated 
w  a  a 

at  the  spatial  frequency  k  •  fT.  The  result  of  applying  a  particular  shading  a^  is  to 

convolve  the  true  (infinite  aperture)  spatial  spectrum  with  the  Fourier  transform  of 

a  ,  m*0, . . . ,  M“l,  thereby  smearing  the  true  spectrum  and  limiting  the  resolution  that  can 

a 

be  obtained. 

HI.  Hlahrtaaolntion  luiglmi 

One  can  obtain  a  set  of  weights  (or,  equivalently,  a  steering  vector)  to  achieve 
better  resolution  by  adapting  them  to  the  particular  noise  field  and  signal  field 
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iapiaiing  on  the  err  my,  la  these  adaptive  beamforming  schemes,  the  steering  vector  is  the 
solution  to  an  optimization  problem[4,  11]:  find  the  steering  vector  which  minimizes  the 
energy  in  the  beam  subject  to  the  constraint  ±'Z  *  1*  where  2  1*  the  constraint  vector. 

The  energy  contained  in  a  beam  can  be  expressed  by  the  quadratic  form  BA  where  £  denotes 
the  empirical  spatial  correlation  matrix  of  the  Fourier  transforms  of  the  sensor  outputs. 
The  correlation  matrix  £  is  usually  estimated  from  the  sensor  outputs  by  a  variation  of 
the  Bartlett  procedure.  The  output  of  each  sensor  is  sectioned  and  windowed.  The  Fourier 
transfozm  of  each  section  is  evaluated  and  the  vector  J^(f^)  of  Fourier  transform  values 

at  the  frequency  fQ  across  the  array  for  the  section  is  formed.  Assuming  that  K 

sections  are  available,  £  is  computed  according  to 

1  •  <s> 
1*0 


Usually  C  is  taken  to  be  the  number  of  statistically  independent  terms  used  in  the 
empirical  computation  of  the  correlation  matrix  £;  K  is  frequently  referred  to  as  the 
time-bandwidth  product. 

The  solution  to  this  optimization  problem  is 


A 


-i  • 

2'*  Z 


(6) 


III.  resulting  nine  of  the  energy  in  the  bean  is 


A'U  “  (Z'*”1!)"’1. 


(7) 


In  the  so-called  high-re solution [5]  or  minimum  energy  scheme,  the  constraint  vector  2  1*  a 

plane-wave  direction-of-look  vector  I,  each  element  of  which  is  given  by  W  *  gj^nmk 

m 

where  k  corresponds  to  a  specific  spatial  frequency.  As  k  *  fT  and  f  is  known,  each  value 
of  k  corresponds  to  a  specific  per-channel  delay.  Defining  0  to  the  bearing  of  the  source 
relative  to  array-broadside#  T  *  (d/fX)sin  C.  Consequently,  spatial  frequency  k  is 
related  to  bearing  0  as  k  *  (d/X)sin  0.  The  constraint  A'X  ■  1  fixes  the  gain  of  the 

steering  vector  in  the  direction-of-look  1  to  be  unity.  Forcing  the  energy  to  be  minimum 
thereby  reduces  the  contributions  from  plane-waves  arriving  from  other  directions  and  from 
the  noise  field.  The  energy  in  the  beam  corresponding  to  the  spatial  frequency  k  is 
expressed  by 


SME(k)  *  d’*”1!)”1. 


(8) 


The  bearing  of  the  target(a)  i»  determined  by  finding  the  spatial  frequency(s)  at  which 
the  quantity  in  equation  (8)  achieves  aaxiaa.  The  naxiaua  likelihood  spectral  estiaate  is 
closely  related  to  the  ainiana  energy  estiaate,  but  differs  fron  it  in  an  iaportant  way. 
The  asxiaua  likelihood  uses  the  noise-only  correlation  aatrix  in  its  evaluation  of  the 

beaa  energy. 
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s„  (k)  -  d'l-1!)-1 

ML  n 


(9) 


Here  R  »  ErNN*  1 .  the  theoretical  correlation  natrix  of  the  noise  component  of  Note 

n 

specifically  that  the  natrix  is  assumed  to  ba  a  known  quantity  and  is  not  computed  from 

empirical  data.  The  maximum  likelihood  method  yields  the  optimum  array  gain[6,  7].  The 
maximum  likelihood  and  minimum  energy  methods  yield  the  same  array  gain  when  the  beam  is 
steered  toward  the  acoustic  sourcetd]. 

H*  l mwiu  Saafllaiiaa 

From  the  viewpoint  presented  here,  there  is  great  flexibility  in  choosing  the 
constraint  vector  £.  One  wonders  if  there  is  a  particular  choice  for  the  constraint 
vector  which  can  maximize  the  spatial  resolution  of  source  bearing.  For  example,  choose  a 
constraint  vector  of  the  form 


Z  -  CW 


(10) 


where  C  is  a  matrix  to  be  described*  The  energy  in  the  beam  when  steered  toward  the 
source  would  then  be 


s(k>  -  <rc*"l<arx.  (11) 

Suppose  the  vector  1  corresponded  to  an  actual  plane-wave  source.  If  C  were  a  matrix 
having  the  property  that  this  choice  of  W  lay  in  the  null  space  of  the  matrix  (Cl  *  0 ) , 

the  energy  in  the  beam  when  steered  in  this  direction  would  be  infinite.  If  the  direction 
vector  corresponding  to  the  plane-wave  source  were  the  only  direction  vector  that  lay  in 
the  null  space  of  C,  one  would  therefore  obtain  a  marked  indication  of  the  bearing  of  the 
source. 

While  it  is  theoretically  possible  to  have  a  perfect  indication  of  source  bearing  by 
this  approach,  the  difficulty  lies  in  finding  the  matrix  C.  This  matrix  must  have  the 
property  that  direction  vectors  lying  in  its  null  space  correspond  only  to  plane— waves 
eminating  from  actual  sources.  Assuming  that  the  bearing  of  the  sources  is  not  known, 
construction  of  the  matrix  C  would  seem  impossible.  However,  one  can  construct  a  matrix 
having  a  null  space  consisting  of  vectors  which  closely  resemble  direction  vectors  of 
source  plane— waves.  The  key  idea  of  this  procedure  is  to  analyze  carefully  the 
eigenvectors  and  eigenvalues  of  the  correlation  matrix. 

IY-A*  UA  Eig.nY.ctor  Method 

The  eigenvectors  of  the  matrix  R  are  defined  by  the  property 

i  “  1 . *  (12) 

where  is  the  eigenvalue  associated  with  the  eigenvector  V^»  As  correlation  matrices 

are  conjugate-symmetric  (Hermitian),  the  eigenvectors  form  an  orthonormal  set.  Assume 
that  the  sensor  outputs  contain  one  signal  and  noise  uncorrelated  with  the  signal 
(X  *  StM)  •  The  result  of  computing  1  according  to  equation  (5)  for  sufficiently  large  I 
is 
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E  -  <r2Q  +  a25g' . 
&  s 


(13) 


Q  is  the  noise  oorrslstion  matrix  normalized  to  have  a  trace  equal  to  the  dimension  M  of 
the  matrix  1*  £  is  the  direetion  vector  of  a  plane-wave  source  and  has  a  squared-norm 

equal  to  M.  The  cross-terms  involving  signal  and  noise  are  assumed  to  be  negligible*  If 
the  noise  correlation  matrix  equals  the  identity  matrix  (i.e.*  only  sensor  noise  is 
assumed  to  be  present)*  the  eigenvector  corresponding  to  the  largest  eigenvalue  (hereby 
referred  to  as  the  "largest  eigenvector")  is  equal  to  the  vector  £  with  eigenvalue  equal 
2  2 

to  a^+Uo^*  The  remaining  M-l  eigenvectors  of  R  consist  of  those  vectors  orthogonal  to  £ 

2 

and  each  has  eigenvalue  or  *  If  p  incoherent*  linearly— independent  signals 

n 

Sj,  i  »  l,...,p,  are  present  so  that  the  correlation  matrix  is  of  the  form 


K 


*2I  + 
n 


1 


,iMl* 


(14) 


the  p  largest  eigenvectors  correspond  to  the  signal  terms  and  the  M-p  smallest  are 
orthogonal  to  all  of  the  signal  direction  vectors.  Mote  that  the  largest  eigenvectors  are 
not  necessarily  equal  to  the  signal  direction  vectors  in  this  case;  these  eigenvectors 
comprise  as  orthonormal  basis  for  the  vector  space  containing  the  signal  vectors* 
Consequently*  one  cannot  always  inspect  the  largest  eigenvectors  and  determine  the  signal 
vectors  directly. 

• 

Define  C_  to  be  the  sum  of  the  outer  products  of  the  M-p  smallest  eigenvectors 


°EV 


(15) 


As  the  p  largest  eigenvectors  are  orthogonal  to  each  of  the  M-p  smallest* 

i  *  M-p+1# . .  .,M.  As  the  p  largest  eigenvectors  span  the  space  containing  the  signal 

vectors*  each  signal  lies  in  the  null  space  of  and  *  0.  In  this  manner* 

perfect  resolution  of  the  bearing  of  multiple  sources  can  be  obtained  from  an  eigenvector 
analysis  of  the  correlation  matrix  R. 

Note  that  in  computing  the  beam  energy  (equation  ID*  the  matrix  need  never  be 
computed.  The  correlation  matrix  R  can  be  expressed  in  terms  of  its  eigenvectors  as 


R 


WL 


(16) 
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(17) 


Because  of  tlxe  orthogonality  property  of  the  eigenvectors  of  a  correlation  matrix,  the 
quantity  C^l  appearing  in  equation  (11)  becomes 


c£v*'1cev 


J  ±L.T. . 

.L,  X .  11 

i»l  i 


(18) 


It  is  this  matrix  which  is  computed  in  the  evaluation  of  the  quadratic  form  of  equation 

(11). 


SEVlk)  ■  <I'CyflcEVirl 

The  following  sequence  of  computations  constitute  the  eigenvector  method. 

1.  Compute  the  correlation  matrix  X. 

2.  Decompose  the  matrix  X  into  its  eigenvectors  and  eigenvalues. 

3.  Determine  the  number  p  of  sources  present  in  the  acoustic  field. 

4.  Compute  the  clergy  in  the  beams  corresponding  to  all  possible  bearings  (equation 
19). 

5.  The  major  peaks  in  this  spectrum  correspond  to  acoustic  sources. 

The  methods  of  0wsley[12],  Bienvenu[2,  3],  and  Schmidt [13]  are  similar,  but  differ 

somewhat  from  the  eigenvector  method  just  described.  The  expression  for  X  ^  is  truncated 
as  in  equation  (5);  however,  the  small  eigenvalues  are  set  to  the  same  value  (taken  here 
to  be  unity).  Instead  of  evaluating  the  quadratic  form  of  equation  (11),  the  spectral 
estimate  of  the  MUSIC  method  can  be  expressed  by 

Wk)  -  (*'CE»CEV*>'1  • 

This  estimate  can  also  ba  evaluated  ia  the  same  maoaer  as  equation  (11)  if  the  matrix  C  is 
radafiaad  to  be 

Smsic  *  1'  VS  til) 

A  steering  vector  can  be  defined  (equation  6)  when  the  latt  *r  approach  is  used  to  express 
the  MUSIC  spectral  estimate.  Under  the  conditions  just  described,  this  spectral  estimate 


IX  ..,.f 
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also  lias  the  capability  of  resolving  source  bearing  perfectly. 

The  preemption  of  the  preceding  analysis  has  been  that  the  noise  correlation  matrix 
Q  is  equal  to  the  identity  matrix.  It  is  this  key  assumption  which  leads  to  perfect 
bearing  resolution.  Imperfection  in  bearing  resolution  occurs  when  this  presumption  is 
false.  The  noise  field  may  contain  more  than  just  sensor  noise:  for  example,  isotropic 
noise  may  also  be  present.  In  addition,  a  finite  amount  of  averaging  is  used  to  compute 
the  correlation  matrix  R.  Even  if  the  noise  field  were  spatially  white,  an  empirical 
noise  correlation  matrix  would  not  be  an  identity  matrix.  Either  of  these  situations  can 
reduce  the  resolution  of  the  eigenvector  and  MUSIC  methods. 

The  performances  of  the  eigenvector  method  and  of  the  MUSIC  method  under  these  more 
realistic  conditions  are  analyzed  mathematically  in  the  appendix.  From  this  analysis,  an 
approximate  lower  limit  on  the  energy  in  the  beam  when  steered  on- target  using  the 
eigenvector  method  is  found  to  be 


S(k) 


(21) 


A  similar  result  is  obtained  for  the  MUSIC  method.  The  quantity  y  is  a  measure  of  the 
spread  of  the  eigenvalues  cu  of  Q.  For  the  eigenvector  method,  this  quantity  is  given  by 


and  for  the  MUSIC  method  by 


2 

ymdsic 


(23) 


2  2 

^ MUSIC  ten**  t0  smaller  than  Ygy*  In  either  case,  the  quantity  y  is  an  implicit 

function  of  the  time-bandwidth  product  E.  Generally  speaking,  y  will  decrease  with 
increasing  I,  tending  toward  the  spread  of  the  eigenvalues  of  the  spatial  correlation 
matrix  of  the  underlying  noise  process.  If  the  theoretical  noise  correlation  matrix 
equals  the  identity  matrix  (i.e.,  spatially  white  noise),  the  spread  of  its  eigenvalues  is 
zero.  Otherwise,  the  spread  is  nonzero. 


The  first  term  in  the  denominator  of  equation  (21)  determines  how  large  the  energy 
peak  will  be  if  an  infinite  time-bandwidth  product  were  available.  Ihe  second  term 
describes  how  the  size  of  the  energy  peak  depends  on  E.  The  larger  of  these  will  dominate 
the  expression  in  equation  (21).  When  steered  off-target,  the  energy  in  the  beam  produced 
by  the  eigenvector  method  coincides  with  that  produced  with  the  minimum  energy  approach. 
The  ratio  (/  of  the  on-target  to  off-target  beam  energy  can  be  used  to  assess  the  size  of 
the  peak  in  the  beam  energy  as  the  direction  vector  1  is  scanned  through  all  possible 
bearings.  When  only  sensor  noise  is  present,  this  quantity  is  given  in  the  minimum  energy 
method  by [10] 
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In  the  eigenvector  aatiiod  and  the  MUSIC  aathod.  this  quantity  is  given  by 


(/ 


(25) 


From  simulation  studies,  this  quantity  can  be  somewhat  larger  (a  few  dB)  in  the  MUSIC 
nethod  than  in  the  eigenvector  aethod.  Conparing  equations  (24)  and  (25),  one  concludes 
tnat  when  the  eigenvector  aethod  is  used,  the  array  appears  to  consist  of  a  number  of 

elements  I  given  by 


X  - 
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In  aost  instances,  this  quantity  is  larger  than  the  actual  number  of  sensors  M. 

IZ.fi.  Simulation  Results 

Computer  simulations  were  used  to  evaluate  the  eigenvector  and  MUSIC  methods  and  to 
compare  them  with  the  ainimum  energy  method.  In  these  studies,  sequences  of  the  data 
vector  fij  •  fi  +  *«re  produced.  The  parameters  of  the  signal  vector  £  were  defined  as 

described  in  section  II.  The  noise  vector  fi,  consisted  of  identically-distributed  complex 

Gaussian  noise  components.  The  covariance  aatrix  of  this  random  vector  could  be  specified 
by  the  user.  Each  noise  vector  was  generated  to  be  statistically  independent  of  all  other 
noise  vectors.  The  empirical  correlation  matrix  R  was  then  computed  as  in  equation  (5) 
and  its  eigenvectors  and  eigenvalues  computed  according  to  a  QL  algorithm [8] . 

With  one  exception,  each  step  of  the  procedure  outlined  above  for  the  eigenvector 
method  was  followed.  As  the  number  p  of  sources  present  in  the  acoustic  field  was  known 
by  the  user,  it  was  supplied  by  him.  In  a  physical  situation,  this  parameter  would  not  be 
so  readily  known.  However,  the  purpose  of  the  simulations  was  to  determine  the  validity 
of  the  theory  and  to  test  how  well  the  methods  could  perform.  The  effects  of  inaccurate 
choices  for  p  are  described  in  a  later  section. 

A  comparison  of  the  beam  energies  produced  by  the  minimum  energy  method  and  by  the 
eigenvector  method  is  shown  in  figure  la  when  one  source  is  present  in  the  acoustic  field. 
A  similar  comparison  is  found  in  figure  lb  for  the  MUSIC  method.  Note  that  the  height  of 
the  main  peak  relative  to  the  background  noise  level  varies  with  E,  the  time-bandwidth 
product,  in  both  methods  whereas  it  does  not  in  the  ainimum  energy  aethod.  In  this 
example,  sensor  noise  (i.e.,  spatially  white)  is  present.  For  the  array  length  used,  the 
second  term  in  denominator  of  equation  (21)  was  larger  than  the  first.  Simulated  and 
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theoretical  value*  of  C^y  are  compared  in  Table  1.  The  theoretical  prediction  of  the 

value  of  this  quantity  is  close  to  that  obtained  from  the  simulat ions. 

The  spectra  obtained  from  the  eigenvector  method  and  the  MUSIC  method  differ  only 
slightly  in  these  examples*  The  latter  method  tends  to  produce  much  flatter  off-target 
spectra*  The  snail  eigenvalues  of  the  correlation  matrix  R  tend  to  represent  the  noise 
field  (see  the  discussion  following  equation  13)  and  the  equalization  of  these  values  vill 
"whiten"  the  background  noise,  thereby  resulting  in  a  flat  spectrun.  This  effect  is 
further  illustrated  by  considering  a  case  where  isotropic  noise  is  present  in  the  acoustic 
field  (figure  2)* 

£[•£'  >*»olotion  si.  Multiple  Targets 

Cox[6]  derive*  analytic  expression*  for  the  limits  to  which  the  minimum  energy  method 
applied  to  a  linear  array  can  resolve  two  equi-strength,  incoherent  acoustic  sources. 
There,  "resolution"  is  defined  as  the  minimum  bearing  separation  at  broadside  at  which  two 
targets  can  be  distinguished  by  evaluating  beam  energy.  Here,  a  slight  dip  in  the  beam 
energy  is  required  when  the  array  is  steered  between  the  sources.  The  critical  factors 
determining  the  resolution  of  beams  formed  by  the  minimum  energy  approach  are  aperture 

(defined  as  the  spatial  extent  of  the  array  relative  to  a  wavelength  -  *  the  number  of 

2  2 

elements  in  the  array  (X),  and  the  sensor  signal-to-noise  ratio  (a^/a^).  Cox's  results 
are  summarized  in  figure  5  of  his  paper[6j;  an  approximation  to  those  results  is 


(27) 


where  0  is  bearing  $ep^rMtion  of  the  two  targets. 

When  the  eigenvector  method  is  spplied  in  situations  such  ss  these,  targets  are  more 
easily  resolved  end  furthermore,  the  resolution  capabilities  of  the  array  are  increased. 
Figures  3  and  4  display  typical  examples  of  these  cases.  A  theoretical  prediction  of  the 
degree  to  which  resolution  is  increased  can  be  obtained  from  equation  (27).  If  one 

substitutes  X  evaluated  by  equation  (2$)  for  M,  the  value  of  0  thus  obtained  is  the 
resolution  limit  of  the  eigenvector  method  and  the  MUSIC  method.  A  comparison  of  the 
resolution  obtained  from  some  of  the  simulations  with  that  predicted  by  the  theory  is 
shown  in  Table  II.  The  degree  of  agreement  between  theory  and  simulation  results  implied 
by  Table  II  is  valid  for  all  of  the  simulations. 

IY.D.  Effect  qf  an  Improper  Choice  for  p, 

The  theoretical  and  simulation  results  presented  thus  far  are  valid  only  when  the 
parameter  p  equals  the  actual  number  of  acoustic  sources.  In  practice,  this  quantity  may 
not  be  known  and  one  questiona  the  sensitivity  of  the  eigenvector  and  MUSIC  methods  to  in 
incorrect  choice  of  p.  This  sensitivity  was  studied  through  the  simulations;  no  analytic 
results  were  obtained  on  this  issue. 

An  incorrect  choice  for  p  has  different  effects  on  the  eigenvector  and  MUSIC  methods. 
In  both  methods,  choosing  p  too  small  does  ppi  result  in  s  beamformer  having  superior 
bearing  resolution  properties  to  the  minimum  energy  method.  In  the  eigenvector  method, 
the  spectra  tend  to  resemble  those  obtained  with  the  minimum  energy  method.  In 

particular,  by  assuming  no  sources  are  present  (setting  p  •  0  or  the  matrix  •  I),  the 

EV 


Johnson  Mmf 


IapmiAg  Bearing  Resolution 


eigenvector  aethod  is  ozmotly  tiie  minima  energy  method.  In  contrast,  the  MUSIC  aethod 
tends  to  produce  a  number  of  spectral  peaks  equal  to  p.  For  example,  if  a  value  of  zero 
is  chosen  for  p,  a  uniformly  flat  spectrum  results.  The  peaks  that  result  from  nonzero 
choices  tend  to  correspond  to  the  bearings  of  acoustic  sources;  however,  which  sources  are 
thus  located  is  not  easily  predicted.  When  p  is  chosen  too  large,  the  bearing  resolving 
capabilities  of  either  aethod  are  not  greatly  reduced.  The  spectra  produced  by  the 
eigenvector  aethod  tend  not  to  vary  from  that  obtained  with  the  proper  value  of  p.  The 
MUSIC  aethod  tends  Co  produce  spurious  peaks  that  do  not  correspond  to  physical  sources. 
The  effects  are  illustrated  in  figure  5. 

V.  Re solution  and  Array  Gain 

While  these  approaches  increase  the  resolution  of  bearing,  this  increase  in 
resolution  is  accompanied  by  a  decrease  in  the  array  gain.  To  show  this,  assume  the 
correlation  aatrix  R  is  of  the  form  given  in  equation  (13).  The  signal* to-noise  ratio  at 

2  2 

each  sensor  is  therefore  o signal-to-noise  ratio  in  the  beam  output  is  the 
quantity 

A'SS'A 

-  (28) 

Hi.  array  gain  6  is  the  ratio  of  these  signal-to-noi.se  ratios. 

A'SS'A 

G  •  — — —  (29) 

A 

Recalling  that  the  steering  vector  in  this  case  is  given  by 

R"1C| 

a  -  - zr-  (30) 
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and  setting  the  direction  vector  to  correspond  to  the  source  (W 
bee ones 


■  1) .  the  array  gain 


A a  the  natrir  Q  is  given  by 


S'  c*  ■T1qr”1C5 


the  denoainator  of  equation  (31)  beooaes 
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so  that  we  obtain 
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2 

-JL 


i'C'£_ICS 


2 

cr 


$ 


For  the  eigeovetor  method. 


V1  ■  ‘'‘Sr  *  CEV*'1,=EV 
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(35) 


and  the  array  gain  be cones 
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(36) 


In  the  MUSIC  method,  the  ratio  appearing  in  the  denominator  of  equation  (34)  can  be 
bounded  using  the  Schwarz  inequality* 

S'^sic^Srosic*  2 

S,<Wsic*'lSS,*‘l|irasicS 

Equality  occurs  only  when  the  M-l  smallest  eigenvalues  of  R  are  identical  (i.e.f  Q  *  I), 

This  bound  can  be  used  in  equation  (34)  to  obtain  an  upper  bound  on  G1ITTC,rri  if  the  bound  is 

MUSIC 

2 

not  saaller  than  a  .  As  the  expression  thns  obtained  equals  G_,  (equation  36),  this 

s  nv 

condition  is  satisfied.  Consequently, 


-L 


(37) 
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Considering  equation  (36),  is  a  aonotonically  increasing  /auction  of  the  quadrstic 

fora  £'Cgp£  Therefore,  whenever  one  decreases  this  quadrstic  fora  to  iaprove  the 

indication  of  bearing  (equation  19),  the  array  gain  decreases  in  the  eigenvector  aethod. 
Because  of  the  relationship  given  in  equation  (38),  the  array  gain  obtained  with  the  MUSIC 
aethod  also  decreases.  In  the  liait,  perfect  indication  of  bearing  (a  zero— valued 
quadratic  fora)  corresponds  to  zero  array  gain  with  either  aethod. 
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II.  Conclosions 

The  eigenvector  method  can  enhance  the  bearing-resolving  capabilities  of  an  array. 
Here,  the  eigenvectors  and  eigenvalues  of  the  correlation  matrix  must  be  found  and  the 
weighted  sua  of  Fourier  transforms  of  the  eigenvectors  computed*  In  the  minimum  energy 
method,  the  inverse  of  the  correlation  matrix  must  be  found  and  the  quadratic  fora  of 
equation  (8)  computed*  Roughly  speaking,  the  computational  complexities  involved  in  the 
use  of  the  eigenvector  method  are  not  excessive  when  compared  to  those  required  in  the 
minimum  energy  method* 

The  degree  to  which  the  resolving  power  is  increased  is  related  to  the  quantity  9* 

1/4 

Because  of  equation  (27) ,  this  increase  is  proportional  to  (H/M)  *  Consequently,  to 

increase  the  bearing  resolution  by  a  factor  of  2  requires  the  virtual  number  of  sensors  9 
to  be  16  times  the  actual  number*  Referring  to  Table  I,  such  large  virtual  array  lengths 
can  be  obtained  only  when  large  time-bandwidth  products  are  possible.  Under  these 
circumstances,  enhanced  bearing  resolution  is  possible.  For  a  given  t im e-bandwidth 
product,  the  smaller  the  number  of  elements  in  the  array,  the  greater  the  increase  in 
bearing  resolution. 

The  eigenvector  method  and  the  MUSIC  method  produce  quite  similar  results.  They 
differ  in  at  least  two  respects,  however.  The  first  is  that  a  non—  zero  value  of  p,  the 
number  of  assumed  sources  in  the  acoustic  field,  must  be  chosen  in  the  MUSIC  method.  If 
the  value  of  p  is  not  close  to  the  actual  value,  the  spectra  thus  obtained  can  differ  from 
that  obtained  with  the  proper  value:  spurious  peaks  appear  and/or  peaks  can  be  missed.  In 
contrast,  the  eigenvector  method  is  less  sensitive  to  the  choice  of  p.  Second,  the  shape 
of  the  spectrum  of  the  background  "noise"  is  drastically  altered  in  the  MUSIC  method.  For 
example,  the  variations  due  to  low-level  sources  or  to  the  physical  noise  spectrum  are 
lost  (see  figure  2,  for  example).  This  portion  of  the  spectrum  can  also  vary  as  p  is 
changed;  this  effect  is  much  less  pronounced  in  the  eigenvector  method* 

The  increase  in  bearing  resolution  is  obtained  at  the  expense  of  array  gain* 
Consequently,  if  more  than  bearing  information  is  required,  other  techniques  should 
probably  be  used  to  obtain  them*  One  can  conceive  of  the  eigenvector  method  being  used  to 
acquire  source  bearing  and  this  information  being  used  to  steer  a  beam  with  the  minimum 
energy  method  so  as  to  analyse  the  waveform  produced  by  the  source*  Mote  that  this  two- 
step  procedure  need  only  be  sequential  in  a  conceptual  manner.  Because  of  the  close 
relationship  between  the  two  methods,  obtaining  the  steering  vector  for  the  minimum  energy 
beamformer  means  including  more  terms  in  the  eigenvector  decomposition  of  R. 

The  decrease  of  array  gain  with  increasing  resolution  raises  many  theoretical  issues. 
The  minimum  energy  method  is  known  to  yield  the  optimal  value  of  array  gain. 
Consequently,  any  method  which  has  greater  resolution  capabilities  cannot  also  increase 
array  gain*  Can  array  gain  be  maintained  while  increasing  resolution  or  is  increased 
resolution  always  obtained  at  the  expense  of  array  gain?  The  present  method  has  the 
latter  property.  A  theoretical  understanding  of  the  limits  to  which  array  gain  and 
resolution  can  be  traded  against  each  other  would  be  of  interest* 

The  main  issue  not  addressed  in  this  study  is  the  determination  of  the  aumber  of 
sources  -  p*  From  the  analysis  presented  in  section  IV,  the  number  of  sources  corresponds 
to  the  number  of  dominant  eigenvalues  in  the  matrix  R.  Determining  p  in  this  way  can  be 
difficult  when  small  time-bandwidth  products  are  available  and  isotropic  noise  is  present. 
Reasonably  accurate  methods  of  determining  p  from  the  eigenvalues  of  R  are  not  known  at 
this  time. 
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Let  tie  vector  J  of  sensor  outputs  be  of  tie  form  l-c^4  +  where  S  denotes  tie 

source  direction  vector  as  before  and  fi  denotes  additive  noise.  Tie  correlation  matrix  R 
is  computed  empirically  according  to  equation  (5)  to  yield. 
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statistical  estimate  of  tie  noise  correlation  matrix  and 


an  estimate  of  tie  average  noise  component.  Tie  vectors  are  assumed  to 


be  statistically  independent  random  vectors;  each  component  of  has  zero  mean  and  unity 

variance.  Consequently*  tie  components  of  tie  vector  3  iave  zero  mean  and  energy  l/I. 
Define  tie  matrix  P  to  be  tie  noise-related  terms  in  equation  (A.l). 


<r  c 

q  +  -*SS'  +  -^5S' 


(A. 2) 


Coaseqncntly,  the  expression  for  K  ia  equation  (A.l)  csa  b.  written  sore  simply  ts 


K 


o2?  +  a 
a 


2 

s 


3£' 


(A. 3) 


The  estimate  of  the  energy  in  the  b.sa  pointed  ia  the  W  dir.ctioa  is  given  by  equation 
(11).  Following  Cox[6],  this  expression  can  be  written  as 


S(k) 


I'C'P^CI 
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JtlMAI 


1+(J)MAXsin2(cI,  I;  P*1) 


(A. 4) 


S  2  2 

where  (”)„._  ■  4'P  &(«/*)  is  the  signal  -to-no ise  ratio  of  tie  beam  output  obtained 
N  MAX.  s  n 

2  -1 

witi  tie  optimally-ciosen  steering  vector  and  sin  (Cl*  4;  P  )  is  tie  sine-squared  of  tie 

angle  between  tie  vectors  Q[  and  4  with  respect  to  tie  matrix  P  lie  matrix  is 

given  by  equation  (15)  wien  p,  tie  number  of  signals  in  tie  acoustic  field,  is  assumed  to 
be  one.  Tie  critical  aspect  of  equation  (A.4)  is  tie  quantity  CgyJ.  This  quantity  can  be 

viewed  as  tie  projection  of  tie  vector  1  onto  tie  set  of  eigenvectors  ortiogonal  to  tie 
largest  eigenvector  of  R«  As  indicated  earlier,  tiese  eigenvectors  are  approximately 

is  orthogonal 

Expression  (A.4) 

therefore  beeoaes 


orthogonal  to  the  signal  weetor.  To  a  good  approximation,  the  rector 

2  -1 

to  1,  thereby  implying  for  all  1  that  sin  (C-~l,  i;  P  )  *  1. 
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When  £  dots  not  correspond  to  the  signal  direction  vector  £,  the  matrix  Cgy  has  little 
effect  on  the  vector  £.  In  this  case*  one  obtains 

2 

c 

S(k)  -  - 

f»P  AI 

the  result  obtained  with  the  minimum  energy  estimate.  Consequently*  one  should  expect  the 
eigenvector  procedure  and  the  minimum  energy  procedure  to  yield  the  same  numerical  values 
when  steered  off* target.  As  the  beam  is  steered  toward  the  source*  the  two  estimates  will 
begin  to  differ  as  the  matrix  C^y  begins  to  sffect  the  vector  V. 

In  the  succeeding  analysis  the  inverse  of  the  matrix  P  is  assumed  to  be  approximately 

equal  to  the  inverse  of  Q  in  the  computation  of  the  quadratic  form  appearing  in  equation 

(A.5).  Inspecting  equation  (A. 2)*  this  approximation  will  be  less  accurate  as  the 

signal-to-noise  ratio  (o  /a  )  increases  and  as  the  amount  of  averaging  (K)  decreases.  The 

s  n 

energy  estimate  can  be  written  approximately  as 
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(A. 6) 


The  quadratic  form  in  equation  (A.tf)  can  be  interpreted  as  the  squared-length  of  the 
vector  Cgyl  with  respect  to  the  norm  induced  by  the  matrix  Q  *.  We  therefore  seek  an 
expression  for  this  quantity  when  £■£. 

Assume  that  the  largest  eigenvector  of  R*  is  given  by  •  £+&  where  ±  is  a 

vector  orthogonal  to  £•  Consider  the  vector  diagram  shown  in  figure  Al.  The  vector  &  is 
defined  to  be  orthogonal  to  the  eigenvector  £+&•  What*  is  sought  in  equation  (A. 6)  is  the 
square  of  the  length  L  of  the  vector  £  projected  onto  As  shown  below*  the  length  of  & 
is  small  compared  to  £;  therefore*  the  quantity  L  will  be  approximately  equal  to  the 
length  of  ±.  To  a  good  approximation*  the  energy  in  the  beam  steered  on-target  is  given 
by 


S(k) 


(A. 7) 


2  —1 
lll.il  .  denotes  tke  square d~norm  of  with  respect  to  <1  . 
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Q  1 


The  definition  of  an  eigenvector  implies  that  this  vector  aaet  satisfy 

ft(fi+£)  -  ^(S+a.)  . 


(A.  8) 


As  the  vectors  £  and  £  are  assnmed  to  be  orthogonal,  the  eigenvalue  1^  can  be  fonnd 
through  the  relationship 

£'ft(£+ft)  -  ^’(i+s.)  -  JUjj 


Using  equation  (A. 3)  for  ft,  we  have 

\  “  m“  + 


ffs®a(£'5tu'£+a,A)  +  m*;. 


(A. 9) 


An  expression  for  the  rector  ±  is  obtained  by  evaluating  the  quantity  After 

some  manipulation  and  assuming  that  the  length  of  £  is  small  compared  to  the  length  of  £, 
we  hare 


t  *  fc.v  -  v1'® 


(A. 10) 


2  -1 

This  expression  for  the  rector  e  consists  of  a  matrix  (X^I-o^Q)  times  the  sum  of 
two  terms*  The  first  term,  denoted  by  is  comprised  only  of  the  signal-related  terms 
and  the  matrix  Q. 

*i  ■  i*-11’ 

Tho  second  tarn  contains  th#  tarns  depending  on  the  average  noise  vector  £. 

£,  ■  Moofl-oo  (£'3)£«  (A.  12) 

l  s  n  s  u 

If  one  assumes  the  noise  vector  3  to  taro  (implying  infinite  statistical  averaging), 

2  -1 

the  vector  £  is  given  by  the  quantity  £  *  (X^I-o^ft)  £^.  Note  that  if  Q  ■  I,  the 

quantity  ■  £  which  implies  £  »  £.  Therefore,  the  signal  vector  £  corresponds  to 

the  largest  eigenvector  of  ft;  this  result  is  consistent  with  the  analysis  described  while 
leading  to  equation  (5).  The  term  £^  oppresses  the  effect  on  the  eigenvector  of  the 
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statistical  averaging  process.  Note  that  if  noise  field  can  be  described  as  containing 
only  sensor  noise  (Q  *  I),  the  expression  for  £  is  dominated  by  fij*  The  vector  can 

be  interpreted  as  the  noise  rector  which  results  when  all  of  its  components  in  the 
direction  of  the  signal  rector  £  are  eliainat ed  (equation  A. 12).  The  squared-magnitude  of 

this  rector  therefore  depends  on  the  "angle"  between  the  rectors  fi  £•  This  angle  will 
be  a  random  quantity  when  the  computation  of  E  is  completed.  This  rector  will  be  largest 

when  5  &nd  &  are  assumed  to  be  orthogonal.  In  this  case,  £'&  *  0,  an  the  expression  for 

fij  becoaes 


*2  "  “Va3* 


(A. 13) 


2  -1 

Define  (S^)  to  b«  the  product  of  (Xgl-o^Q)  and  The  vector  ±  is  therefore 

written  as 


^  "  *1  +  *2 


(A. 14) 


An  expression  for  the  mom  of  jg  can  now  be  obtained.  Let  denote  an  eigenvector  of 
the  matrix  Q  and  the  associated  eigenvalue.  As  these  eigenvectors  are  orthonormal,  any 
rector  can  be  expressed  as  a  linear  combination  of  them.  As  is  also  an  eigenvector  of 
2  -1 

(Xgl-o^Q)  ,  and  c*n  be  written  as 


% 


(°i 


-  i  k 


o.o 

M  in 


(A. 15a) 
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...  3-SS 


Sl~ai°a 


2  °i 


(A. 15b) 
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The  quantity  c*a  be  siaplified.  Using  equation  (A. 9)  and  assuming  both  ±  and  £ 

are  small  compared  to  £,  we  hare 


-  a^)a 


2 

n* 


The  quantities  within  the  parentheses  are  comparable  and  further  assuming 


we 


hare 
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2  2 

X^-a.o  -  Ma  « 
Min-  $ 


(A.ltf) 


The  norm  of  g,  depends  upon  the  angle  between  these  two  vectors*  Based  on  statistical 
arguments*  a  zero-mean  vector  obtained  from  averaging  in  this  case)  will  be  nearly 
orthogonal  to  any  fired  vector*  Consequently , 


llll2  -  IIS,  II2.,  *  IIS,  II2 

Q  Q  « 


(A. 17) 


NS,"2., 


<«i  -  i  Kli-s.|2>: 


* .  2,2 

ttiaM"ai®a) 


To  evaluate  this  expression*  the  relationship  between  the  signal  direction  vector  S  and 
the  eigenvectors  of  4  must  be  specified.  If  £  were  proportional  to  an  eigenvector  of  4* 

the  quantity  IIS1  II  .  would  be  zero.  As  the  norm  of  BL  will  appear  in  the  denominator  of 
1  4  1 

the  expression  for  beam  energy  (equation  A. 7)*  one  can  obtain  a  lower  bound  on  the  energy 
in  the  beam  by  assuming  the  largest-possible  value  for  its  length.  To  approximate  the 

maximal  length  of  assume  that  3  does  not  prefer  anv  of  the  eigenvector  directions  of 

Q.  A  reasonable  mathematical  description  of  this  situation  is  that  m  ai* 

instance*  we  have  using  equation  (A.  16)  that 


"Ml*  -  '-i  i  1  <«i  -  i  K>: 

0  a  Br  i  m 


(A. 13) 


Hi.  quantity  in  tie  summation  d.p.nds  only  on  the  eigenvalues  of  the  matrix  Q.  Define  the 


quantity  tgy  to  be 


4  - 


(A. 19) 


Then,  equation  (A.18)  becomes 


4  2 

n«  m2  .  !a.Iex 

"“l"  -1  4  M  * 

Q 


(A. 20) 
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As  the  ten  2^  is  a  randoa  quantity,  its  squared-nora  is  defined  to  be 


"m;-* 


«:  -  Ifi'Ej 


21 


5 


2  _ 
a  i 
s 


(A. 21 ) 


where  E[*]  denotes  expected  value*  To  a  good  approxinat ion,  is  an  eigenvector  of  the 
correlation  natrix  associated  with  the  randoa  vector  5*  Consequently, 


so  that 


EllE'fiJ2] 


a. 


I 


(A. 22) 


Substituting  equations  (A. 18),  (A. 22).  and  (A. 17)  into  equation  (A. 7),  we  have  finally 


SEV(k) 


2  2 

!»!« I 

2  N 


M 

c 


(A. 23) 


as  an  expression  for  tha  energy  in  the  been  when  steered  toward  the  source. 

The  analysis  for  the  MUSIC  nethod  differs  only  in  detail  frost  that  just  described. 
In  this  nethod,  the  spectral  estimate  is  giwen  by 

Wl>  ’  <*-CEVCEV*>'1  <*•»’ 

Off  target,  Cg^l  approxiaately  equals  5,  implying  that  ^^ujsIC M  ^ea  steered  on 

target,  the  expression  for  the  MUSIC  estiaate  differs  little  froa  that  given  in  equation 
(A«7),  The  significant  difference  is  that  the  norm  of  j:  is  eoaputed  with  respect  to  the 

identity  natrix  instead  of  Q  •  The  quantity  of  interest  is  therefore 

llsll2  -  IIEJI2  ♦  llfijll2* 

Tha  non  of  rsspaot  to  I  equals  that  cos^utsd  with  rsspact  to  Q  th*  sxprsssion 
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foe  the  acta  of  i* 


11^11 


2 


Therefore,  the  on- tar got  beam  energy  in  the  MUSIC  nethod  is  given  by 


SM7SIC(k) 


2  2 

la  lassie  m 

2  M  I 

<J 


where 


MUSIC 


is  defined  to  be 


2 

7JflJSIC 


1 

M 


tti  (ai 


1 
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(A.25) 
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Table  I*  Empirical  and  Theoretical  Peak-to-Background  Ratios 

Hie  results  of  computer  simulations  and  theoretical  predictions  of  the  value  of  (J  are 
shown.  In  the  simulations*  an  equally-spaced  linear  array  containing  10  sensors  was 
assumed  to  be  present  in  an  acoustic  field  containing  sensor  noise  and  a  plane-wave 
source*  The  sensor  signal-to-noise  ratio  was  0  dB.  Values  of  y  were  computed  in  separate 
simulations  from  noise-only  correlation  matrices  having  the  same  time-bandwidth  product. 
The  source  was  assumed  to  be  narrowband*  with  all  of  its  source  energy  concentrated  in  one 

temporal-frequency  analysis  bin.  The  sensor  spacing  is  one-half  wavelength.  Plot  of  beam 

/\ 

energy  vs.  bearing  were  obtained  and  empirical  values  of  (/  estimated.  The  quantities  Um 

ME 

A 

and  (/gy  correspond  to  these  empirical  values. 


A  A 


K 

H 

</me  (dB) 

%  (dB) 

20 

18 

10 

13 

13 

50 

46 

11 

18 

17 

100 

90 

12 

22 

20 

200 

183 

11 

27 

23 

500 

461 

11 

29 

27 

1000 

940 

11 

32 

30 

2000 

1784 

11 

34 

33 

5000 

4597 

11 

39 

37 
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le  II.  Empirical  and  Theoretical  Resolution  Liaita 

ra suits  of  computer  simulations  and  theoretical  predictions  of  resolution  limits  are 
marized.  Hie  array  configuration  used  in  Table  I  vas  used  here.  The  sources  were 
metrically  looated  about  broadside  (0  degrees).  Plots  of  bean  energy  vs.  bearing  were 
aimed  and  the  separation  between  the  sources  reduced  until  they  could  just  be  resolved. 

measurements  of  separation  were  made  in  half-degree  increments.  The  angular 
entities  6  are  indicated  in  degrees. 
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Eina  fimiaai  * 

Figure  i.  Energy  in  beans  formed  by  the  minima  energy  method  (ME)  and  by  the  eigenvector 
method  (EV)  are  plotted  against  bearing  for  a  linear  array.  Energy  is  expressed  in  dB 
relative  to  the  peak  value.  Bearing  is  expressed  in  degrees  with  zero  corresponding  to 
broadside.  The  sensors  are  equally  spaced  and  separated  by  half  a  wavelength.  The  source 
was  assumed  to  be  narrowband,  with  all  of  its  source  energy  concentrated  in  one  temporal— 
frequency  analysis  bin.  In  each  sub-figure,  sensor  noise  and  one  source  located  at  0  are 
present  in  the  sound  field;  the  sensor  signal-to-noise  ratio  is  0  dB.  The  results 
obtained  when  two  time-bandwidth  products  (K)  are  used  are  shown  in  each  sub-figure  for 
each  nethod. 

Figure  £.  Bean  energy  is  plotted  against  bearing  when  one  source  (located  at  0°)  and 
isotropic  noise  are  present  in  the  sound  field.  The  array  conf iguration  is  similar  to 
that  described  in  figure  1,  the  only  exception  being  that  the  sensor  spacing  is  three- 
eighths  of  a  wavelength.  The  sensor  signal-to-noise  ratio  is  0  dB.  The  results  of 
applying  the  minimum  energy  (ME),  eigenvector  (EV).  and  MUSIC  methods  are  shown.  The 
time-bandwidth  product  here  is  50;  the  theoretical  values  of  (J  corresponding  to  this 
situation  are  (/^  *  15.7  dB  and  *  16*4  dB. 


Figure  3..  Beam  energy  is  plotted  against  bearing  when  two  sources  are  present  in  the 
acoustic  field.  The  conventions  defined  in  Figure  1  apply  to  this  plot.  Here,  the  sensor 
si^aal-to-noise  ratio  of  each  source  is  0  dB  and  the  the  sources  are  located  at  -5  and 


Figure  Beam  energy  is  plotted  against  bearing  when  two  sources  are  present  in  the 
acoustic  field.  The  conventions  defined  in  Figure  1  apply  to  this  plot.  Here,  the  sensor 
signal-to-noise  ratio  of  each  source  is  0  dB  and  the  the  sources  are  located  at  -3°  and 


Figure  £.  Beam  energy  is  plotted  against  bearing  with  the  number  p  of  terms  truncated 

from  the  eignevector  expansion  of  R  *  as  a  parameter.  The  parametric  beam  energy 
functions  in  each  panel  is  plotted  with  the  same  vertical  scaling.  A  linear  array  of  ten 
equally-spaced  sensors  (spacing  equal  to  three-eights  of  a  wavelength)  is  present  in  an 
acoustic  field.  Three  incoherent  sources  are  present  in  the  acoustic  field:  two  have 
unity  amplitude  and  are  located  at  bearings  +5°  and  -5°  while  the  third  has  amplitude  of 
one— half  and  bearing  -40  •  Isotropic  noise  is  also  present  in  the  acoustic  field;  the 
sensor  signal-to-noise  ratio  (relative  to  the  larger  signals)  is  0  dB.  The  time-bandwidth 
product  in  both  panels  is  50.  The  upper  panel  displays  the  result  of  applying  the 
eigenvector  method  and  the  bottom  panel  illustrates  the  result  of  applying  the  MUSIC 
method  for  the  same  set  of  data.  Note  that  the  proper  value  of  p  for  these  data  is  p»3. 


Figure  Relationship  of  the  Signal  Vector  and  the  Largest  Eigenvector 
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